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J5 ABSTRACT 

Aims. We investigate the influence of ram-pressure stripping on the internal gas kinematics of simulated spiral galaxies. Additional 
' emphasis is put on the question of how the resulting distortions of the gaseous disc are visible in the rotation curve and/or the full 2D 

CN) I velocity field of galaxies at different redshifts. 

Methods. A Milky- Way type disc galaxy is modelled in combined N-body/hydrodynamic simulations with prescriptions for cooling, 
star formation, stellar feedback, and galactic winds. This model galaxy moves through a constant density and temperature gas, which 
has parameters similar to the intra-cluster medium (ICM). We study five different configurations, in which the direction of the ram 
pressure on the gaseous disc varies from face-on to edge-on. Rotation curves (RCs) and 2D velocity fields of the gas are extracted 
^ , from these simulations in a way that follows the procedure applied to observations of distant, small, and faint galaxies as closely as 

5^ . possible. 

' Results. We find that the appearance of distortions of the gaseous disc due to ram-pressure stripping depends on the direction of 

' the acting ram pressure. In the case of face-on ram pressure, the distortions mainly appear in the outer parts of the galaxy in a very 

symmetric way. This distinction between undisturbed inner part and symmetrically disturbed outer regions is visible in both the 
rotation curve and the 2D velocity field. In contrast, in the case of edge-on ram pressure we find stronger distortions. In particular, a 
mismatch between kinematic centre and the centre of the stellar disc becomes observable. In the RC, there are basically two features 
^ ' visible. The RC is asymmetric as the kinematic centre does not coincide with the optical centre of the galaxy. Secondly, the outer parts 

c 2 ^ , of the RC are declining. The 2D velocity field also shows signatures of the interaction in the inner part of the disc. At angles smaller 

0^ ■ than 45° between the ICM wind direction and the disc, the velocity field asymmetry increases significantly compared to larger angles. 

' Compared to distortions caused by tidal interactions, the effects of ram-pressure stripping on the velocity field are relatively low in all 

CO ' cases and difficult to observe at intermediate redshift in seeing-limited observations. 

CO ' Key words. Galaxies: interactions - intergalactic medium - Galaxies: kinematics and dynamics - Methods: numerical 
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• • 1. Introduction processes, which are partly poorly understood. Another method 
^ ■ , , . , . , , , , , . is to study the total gravitational potential of galaxies and its pos- 
It has been pointed out ever since the late 1970 s that there is an ^.^j^ distortions due to interactions via their internal kinematics, 
evolution of galaxy cluster members with redshift. The fraction Technically it became feasible to observe the full 2D velocity 
^ of star-forming and post-star-forming systems, for example, has ^^j^ ^yp^ j^^^j j^^^^^ -^^^ wavebands using integral 
been proven to increase sipificantly with z (Dressier et al. 1987, ^^j^ ^^.^^ ^jp^^^ ^^^^ SAURON (e.g. Ganda et al. 2006) or 
1999). Consistendy, Butcher & Oemler (1978) reported a strong p^bry-Perot interferometry (e.g. Chemin et al. 2006; Garrido et 
evolution from bluer to redder colours in cluster galaxies, detect- 2002). For intermediate and high redshift galaxies, however, 
ing an excess of blue objects at z=0 5 with respect to lower red- ^j^^^.^ ^^.^ observational studies of 2D velocity fields avail- 
shift systems (this observational finding is usually referred to as ^^j^ ^^^^.^^ ^^^^ ^^^^^.^^ ^^^.^^^ ^j^^.^^ ^^^^ 
Butcher - Oemler effect ). Various physical mechamsms were ^j^^ spectroscopy (e.g. Weiner et al. 2006; Moran et al. 2007; 
proposed which act on the star formation activity and on the mor- ^^^^ ^ ^iegler 2007). To identify distortions and irregularities 
phology of galaxies in clusters besides the hierarchical assembly ^^j^^ ^^^^ ^^^^^ ^^-^-^^^ ^ desirable goal 
of structures (e.g. ram-pressure stripping: Gunn & Gott 1972 - ^^^j^ disentangle different interaction processes due to 
harassment : e.g. Moore et al 1998 - strangulation or star- ^j^^^^ inMence on the velocity field. Numerical simulations offer 
vations : e.g. Larson et al. 1980). We still need to understand the ^ jbiiity to study individually the effects of the various pro- 
specific importance of each of the proposed processes. For that ^^^^^^ Kronberger et al. (2006, 2007) investigated in this way 
purpose It IS crucial to identify and diseritanglediff^erent inter- ^j^^ ^^^^^^ galaxy-galaxy mergers and of tidal interactions 
actions, also at intermediate and high redshift. One possibility ^^^^^^^ j^^^^^ -^^^^^^^ kinematics of model galaxies. 
IS to study the stel ar population of the galaxies which, how- j^^^^ simulations were also used to identify possible observa- 
ever, depends heavily on the current star formation activity. The ^.^^^^ ^.^^^^ observations of the velocity field of distant galax- 
time-dependent star form ation encompasses several non-hnear -^^ (Kapferer et al. 2006; Kronberger et al. 2007). Although de- 

e , jr ■ , , , T u -1 pendencies on the viewing angle and on the spatial resolution 
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troduce nonaxisymmetric and nonbisymmetric features (see also 
Rubin et al. 1999). 

Ram-pressure stripping is another important cluster-specific 
interaction process. The space between the galaxies in a galaxy 
cluster is filled with a hot (~ 10^ K), thin (~ 10^ ions/m^) 
plasma, the so called intra-cluster medium (ICM). This plasma 
exerts a pressure on the inter-stellar medium (ISM), which can 
remove gas from the disc if the force due to ram-pressure strip- 
ping exceeds the restoring gravitational force of the galaxy. This 
mass loss due to ram-pressure stripping was modelled by many 
groups (e.g. Abadi et al. 1999; Vollmer et al. 2001; Roediger & 
Hensler 2005; Quilis et al. 2000; Mori & Burkert 2000; Toniazzo 
& Schindler 2001; Schulz & Struck 2001; and Jachym et al. 
2007). Recently also the effects of the external ram pressure 
on the star formation rate of the stripped galaxy were studied 
(Kapferer et al. 2008, Kronberger et al. 2008). Here we study 
the influence of ram-pressure stripping on the internal gas kine- 
matics of simulated spiral galaxies. 

This interaction mechanism has little influence on the stellar 
structure and could, therefore, be hidden in a galaxy that has a 
regular appearance when only direct imaging is used to search 
for interaction processes. However, it has a strong effect on the 
gas clouds within a galaxy disk disrupting its rotational motion, 
which can be traced by emission lines originating in these gas 
clouds. If spectra of observed galaxies have sufficient signal in 
the continuum, absorption lines can be used to also derive stellar 
rotation curves. Discrepancies between these and the gas motion 
are then strong indications of ram pressure. The velocity field of 
faint distant galaxies is in most cases measured only by emission 
lines (e.g. Ziegler et al. 2002, Bamford et al. 2005, Metevier et 
al. 2006) so that one needs a very thorough understanding of the 
effect of ram-pressure stripping. 

Observable signatures of ram-pressure stripping were also 
studied numerically by Roediger & Briiggen (2006, 2007) using 
Eulerian hydrodynamics. They focussed on the gas distribution 
in the disc and in the wake and thereby tried to deduce informa- 
tion on the interaction geometry. The velocity field in the wake 
has also been investigated. They have, however, not considered 
the velocity field in the disc in detail, as we do it in the present 
work. In a series of papers Vollmer et al. have used sticky particle 
simulations to model specific observed systems (Vollmer et al. 
2000, 2003, 2006, 2008). They have tried to match observed HI 
velocity fields and gas distributions to infer information on the 
interaction scenario. Consequences of ram-pressure stripping for 
observable ICM properties are discussed in Schindler & Diaferio 
(2008) and references therein. This present work goes beyond 
these previous studies by systematically investigating the effects 
of ram-pressure stripping observable with rotation curves from 
optical long-slit spectroscopy and 2D velocity fields at different 
redshifts. 

We also study whether the deviations from a normal rotation 
caused by ram-pressure stripping have a different form than in 
the case of other interactions, so that they could be taken as an 
indicator of what mechanism may be acting on a galaxy. In the 
cluster and group environment it is expected that ram-pressure 
stripping is one of the main drivers for evolutionary effects on 
galaxies. It could be responsible for the global suppression of 
star formation in cluster galaxies (either through an initial star 
burst triggered by gas compression, as found e.g. by Kronberger 
et al. (2008) and Kapferer et al. (2008), or through the subse- 
quent removal of gas from the disk). On the other hand, it can- 
not be responsible for structural changes needed for a complete 
transition from one galaxy type to another. 



Table 1. Particle numbers and mass resolution used for the 
model galaxy. 





particle number 


mass resolution 
[M0/particle] 


Dark matter halo 


300000 


3.5 X 10" 


Disk collisionless 


200000 


1.0 X 10^ 


Gas in disk 


200000 


3.4 X lO'' 


ICM 


1000000 


1.4x 10« 



The paper is organised as follows: in Sect. 2 we present the 
simulation setup used for this work. In Sect. 3 the extraction of 
the observable quantities (i.e. the rotation curves and the 2D ve- 
locity fields) are described. The results are presented in Sect. 
4 and subsequently discussed with respect to previous observa- 
tional and numerical work in Sect. 5. We end with a summary of 
the main conclusions in Sect. 6. 

2. Simulations 

The simulations were carried out with the N-body/SPH code 
GADGET-2 developed by V. Springel (see Springel 2005 for 
details). This code treats the gas of the galaxies and the ICM 
hydrodynamically via smoothed particle hydrodynamics (SPH, 
Gingold & Monaghan 1977; Lucy 1977) while the collisionless 
dynamics of the dark matter and the stellar component is cal- 
culated using an N-body technique. Prescriptions for cooling, 
star formation (SF), stellar feedback, and galactic winds are in- 
cluded as described in Springel & Hernquist (2003). The mass 
loss of the galaxy due to galactic winds, M^,, is assumed to be 
proportional to the star formation rate ^sfr, i-C- M„ - tj^'sfr 
with rj - 2, consistent with the observations of Martin (1999). 
Additionally, the wind contains a fixed fraction x of the super- 
nova energy, which is assumed to be ;t'=0.25, as in Springel & 
Hernquist (2003). 

The initial conditions were built according to Springel et al. 
(2005), based on the analytical work of Mo et al. (1998). The 
model galaxy, which we use, represents a Milky-Way type spi- 
ral galaxy. The total mass of the model galaxy is 1.09x10'^ 
Mq, where the initial total gas mass is 6.8x10^ Mq. The ini- 
tial disc scale length is r,; ~ 3.3 kpc. Compared to previous work 
(Kapferer et al. 2005, 2006; Ki-onberger et al. 2006, 2007), where 
we used a similar sized galaxy, we have significantly increased 
the number of particles, increasing the mass resolution of the 
simulation. Important numerical quantities, such as the parti- 
cle numbers, are summarised in Table [T] For the intra-cluster 
medium we use one million SPH particles, which we distribute 
homogeneously over a volume of 1 Mpc^ with a mass density of 
1 x 10"^** g/cm^ and a constant temperature of 3 keV (~ 3.6 x 10^ 
K). Such a simplified ICM distribution allows the effects of ram- 
pressure stripping to be studied in a clean way, i.e. with as little 
degeneracies of different effects as possible. Additional effects 
from varying density and temperature in the ambient medium 
will be investigated in an upcoming work. 

Starting from these initial conditions we calculate the evolu- 
tion of the model galaxy moving with a constant velocity of 1000 
km/s through the ambient medium for 1 Gyr. We calculated five 
different interaction geometries, where we varied the angle be- 
tween the ICM wind direction and the plane of the disc from 0° 
to 90°. The edge-on and face-on stripping configurations are the 
two extreme cases where the perturbation acts either perpendic- 
ular to the galactic disc or purely in the plane of the gaseous disc. 
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We study these two cases in detail. The other three simulations, 
where this angle is 23°, 45°, and 68°, respectively, are used to 
study how the asymmetry of the velocity field depends on the 
ICM wind direction. 

3. Extraction of realistic rotation curves and 2D 
velocity fields 

In order to be able to compare our results to observations of dis- 
tant, small, and faint galaxies we extract rotation curves and 2D 
velocity fields as described in Kronberger et al. (2006, 2007) in 
a way which closely resembles the observational approach (e.g. 
Ziegler et al. 2003). For the construction of realistic 2D veloc- 
ity fields, we project all gas particles onto a Cartesian, equally 
spaced grid. The spacing is chosen such that the spatial resolu- 
tion corresponds to the angular resolution of current state-of-the- 
art observations. We adopt an angular resolution typical for IFU 
or FPI observations, namely 0.4" (e.g. Chemin et al. 2006). The 
angular resolution of SAURON, for example, would be 0.3" or 
0.9". We calculate the physical resolution according to the given 
angular resolution for a particular redshift using the concordance 
cosmological model. The velocity field of the galaxy is binned 
using this spatial resolution. 

For observations of galaxies at intermediate and high red- 
shift seeing plays a crucial role, as it typically exceeds the an- 
gular resolution of the instrument. To simulate seeing effects on 
the velocity-field measurements we apply a convolution with a 
Gaussian point spread function. A value of 0.8" for the FWHM 
of the Gaussian seeing was adopted, which is a typical value for 
ground based observations (see e.g. Jager et al. 2004) and is rea- 
sonable for the achieved spectral resolution (and a slit width of 
1"). 

From this velocity field we can further extract a rotation 
curve by modelling a slit which is placed over the major axis 
of the galaxy with a typical slit width of 1", as detailed in 
Kronberger et al. (2006). 

4. Results 

In the following we investigate observable projections of the full 
3D velocity field of the simulation, i.e. the rotation curves and 
the 2D velocity fields. Special emphasis is put on the question 
how the given distortions of the velocity field appear at interme- 
diate redshift. 

We quantify the degree of asymmetry of the RC shape, fol- 
lowing an approach used in Dale et al. (2001) and subsequently 
for simulated galaxies in Kronberger et al. (2006). The area be- 
tween the kinematically folded approaching and receding halves 
is divided by the average area under the RC: 

V \mr)\-\V(-r)\\ 
Asymmetry - > 

^ Vo-2(r) + crH-r) 

Here V{r) is the velocity and crir) is the uncertainty of the 
rotational velocity at position r. As in Kronberger et al. (2006) 
we adopt for cr(r) a value of 20 km/s, which is typically the min- 
imum uncertainty in observations of distant galaxies due to the 
limited spectral resolution. 

In order to quantify the distortions and to interpret the 
partly complex structures in the 2D velocity fields we use as in 
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Kronberger et al. (2007) the kinemetry package of Krajnovic et 
al. (2006). This method is based on the assumption that the mean 
velocity along best fitting ellipses is reproduced by a simple co- 
sine law, i.e. 

Via, T) = Vo + VM)cos'i>, (2) 

where a is the length of semi-major axis of the ellipse, ^ for 
discs is the azimuthal angle measured from the major axis in the 
plane of the galaxy. The position angle F and the axial ratio (q - 
bid) of the ellipses are calculated as a function of radius from 
the galactic centre. Deviations from the cosine law are measured 
using an harmonic expansion along the ellipses, i.e. 

N 

y(a, kn{a)cos[nQ¥ - (p„(a))], (3) 

n=l 

where (p is the phase coefficient. Higher order Fourier terms and 
radial changes of F or <7 quantize deviations in the velocity field 
from a simple rotation. Sometimes a simple sine correction for 
the inclination of the disc is applied in observations. As this cor- 
rection simply alters the absolute value of the velocities, we do 
not apply it to our model velocity fields. All VFs were extracted 
at the same inclination i=35°. 

In Fig. [T] we present the 2D velocity field and a long-slit RC 
from an undisturbed model galaxy as seen at redshift z - 0.1. 
For the appearance of the VF at different redshifts we refer to 
Ki-onberger et al. (2007). 

4.1. Rotation curves of ram-pressure stripped galaxies 

First we study the RC of a model galaxy moving/ace-on through 
the ambient medium. After 100 Myr of interaction a clear strip- 
ping radius is visible, which stays roughly constant over time 
and agrees with analytical estimates (Kronberger et al. 2008). 
The stripping radius in ram-pressure affected galaxies is defined, 
as the distance from the galactic centre, outside which the inter- 
stellar medium (ISM) cannot be prevented from being stripped 
by the galactic gravitational potential. We extract rotation curves 
from the galaxy before the clear emergence of a stripping radius, 
i.e. in the early phases of stripping (after ~ 50 Myr), and there- 
after (after ~ 100 Myr). The RC of the isolated model galaxy has 
the typical shape, rising in the inner and flattening in the outer 
parts. After 50 Myr of acting ram pressure first effects of the 
external pressure are visible in the distribution of the gas. The 
strength of the disturbances depends on the local gravitational 
potential. They are strongest at large galactic radii, where the 
restoring force acting against ram pressure is weakest. Whether 
these distortions are visible in the rotation curve depends on the 
radius to which the RC can be measured. We assume that the 
rotational velocity can be determined accurately out to 4 times 
the radial disc scale length of the galaxy. This is a typical value 
obtainable at intermediate redshift (e.g. Ziegler et al. 2003). In 
the top panel of Fig.|2]we present the corresponding RC, which 
does not show strong peculiarities (asymmetry measure ~6%). 

Local measurements extend to much larger radii. We con- 
sider the case where the rotational velocity can be measured out 
to 6 times the radial disc scale length. Already after 50 Myr, the 
gas at such large radii gets a dominant velocity component an- 
tiparallel to the galaxies' direction of movement, i.e. away from 
the galactic disc. How this change of the velocity vector in the 
stripping region translates to a change in the RC of the galaxy is 
shown in the lower panel of Fig.|2l In the outer parts of the RC, a 
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Fig. 1. 2D velocity field and long-slit rotation curve of an undis- 
turbed model galaxy as seen at redshift z - 0.1. The figure of 
the 2D velocity field is taken from Kronberger et al. (2007). 
Overlayed as contours are the best fitting ellipses from the 
kinemetry analysis. 



turn-over of the rotational velocity on both sides to higher nega- 
tive values can be observed. This is a consequence of gas moving 
away from the galactic disc due to the ram-pressure. The asym- 
metry measure of this disturbed RC is ~ 22%, which is low com- 
pared to distortions caused by tidal interactions (see Kronberger 
et al. 2007). Note that the 'direction' of the distortion depends on 
the movement of the galaxy with respect to the observer. If the 
galaxy was moving away from the observer, the rotational ve- 
locities in the outer parts of the RC would go to higher positive 
values. 

After 100 Myr the stripping radius at r~12 kpc is clearly vis- 
ible. Hence also the RC distortion is now visible in this region 
of the disc, and also observable at intermediate redshift. The RC 
for this time step is presented in Fig.|5](top panel). There is a de- 
pendence of the RC shape on the inclination of the galactic disc 
but the important characteristics of the RC remain observable for 
most inclinations. These characteristic features in the outer parts 
of the RC are less prominent if the galaxy is observed almost 
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Fig. 2. Rotation curve of a model galaxy after 50 Myr of ram 
pressure acting face-on. The RC is assumed to be observable out 
to 4 (top panel) and 6 (lower panel) times the radial disc scale 
length, respectively. In the top panel, the RC is similar to the ini- 
tial undisturbed configuration. In the lower panel signatures of 
ram-pressure stripping are clearly visible in the outer parts. Note 
that the 'direction' of the distortion depends on the movement of 
the galaxy with respect to the observer. If the galaxy was mov- 
ing away from the observer, the rotational velocities in the outer 
parts of the RC would go to higher positive values. 



edge-on. There the gas velocity is dominated by a component 
perpendicular to the galactic disc due to the acceleration by the 
ICM wind. If the galaxy is seen edge-on, only the altered disc 
velocity component is observed. The dependence on the viewing 
angle is not as strong as for tidal interactions (see Kronberger et 
al. 2006) as ram-pressure stripping is an axisymmetric distor- 
tion. The characteristic features are observable for all viewing 
angles. 

Next we investigate the RC of a model galaxy flying edge-on 
through the ambient medium. In this case the perturbation acts 
in the plane of the disc. The gaseous disc gets strongly com- 
pressed, which alters the radial pressure gradients and hence the 
rotational velocity of the gas. We find stronger distortions than 
in the face-on case, especially a mismatch between kinematic 
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30 kpc 

Fig. 3. The gas distribution is overlayed as contours on the 
smoothed stellar distribution after 500 Myr of ram pressure act- 
ing edge-on. The white arrow in the top right corner of the figure 
indicates the ICM wind direction. The galaxy rotates counter- 
clockwise. 



centre and the centre of the stellar disc is observable. The off- 
set of the gas kinematic centre is a; 2 kpc after 500 Myr of ram 
pressure acting. In Fig.[3]we show the contours of the gas distri- 
bution overlayed on the smoothed stellar distribution after 500 
Myr of ram-pressure acting (we chose this evolved time step, 
compared to the face-on case, as the distortions are best visible 
there. Before, the distortions are qualitatively the same, but less 
prominent). In contrast to the face-on stripping, there is no dis- 
tinct centre or symmetry visible in the gas distribution anymore. 
In the rotation curve, there are basically two features visible. As 
the kinematic centre does not coincide with the optical centre of 
the galaxy anymore, the RC is slightly more asymmetric than in 
the case of face-on ram pressure (asymmetry measure ~30%). 
Secondly, the RC is declining in the outer part on both sides. 
Note that this feature is in this case only visible in the RC of the 
gas, as only this component is affected by the ram pressure of 
the ambient medium. The collisionless components of the galaxy 
(i.e. the stars and the dark matter) are not or due to the changing 
gas density distribution only mildly affected. If a declining RC is 
also observed in the stellar RC, a tidal interaction affecting also 
the stars and the dark matter is most likely the cause. There is a 
dependence on the viewing angle and the inclination of the disc, 
not strongly affecting, however, the main feature of the disturbed 
RC, namely the declining outer parts. 

For better visual judgement of the asymmetry of the RC, 
both its sides are sometimes plotted on one side. In Fig. |5] we 
present the RCs in such a way for the model galaxy after 100 
Myr of face-on ram pressure and after 500 Myr of edge-on ram 
pressure, respectively. 
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Fig. 4. Rotation curve of a model galaxy after 500 Myr of ram 
pressure acting edge-on. 



In both cases, edge-on and face-on, the long-slit RC can- 
not be used to estimate parameters for a Tully-Fisher study. 
The kinematics of the gas in the galaxy is affected by ram- 
pressure stripping. The unaffected stellar kinematics could be 
used instead to estimate the maximum rotational velocity Vmax- 
However, the luminosity of the galaxy is also affected by ram- 
pressure stripping. In a recent work, we found an enhancement 
of the star-formation rate due to ram pressure (Kronberger et 
al. 2008). A subsequent brightening of the ram-pressure affected 
galaxy can cause an increased scatter in the Tully-Fisher rela- 
tion. 

In Fig.|6]we study how the asymmetry of the VF depends 
on the angle between the ICM wind direction and the plane of 
the disc. The asymmetry measure depends on the viewing angle 
and the inclination, so we averaged over three different lines- 
of-sight and two observed disc inclinations. A clear step in the 
asymmetry measure can be seen at 45°, where the asymmetry 
increases significantly. The scatter, however, increases as well, 
which makes it difficult to distinguish further between different 
ICM wind directions with an observed VF asymmetry. A high 
VF asymmetry is nevertheless a good hint, for a small angle be- 
tween the ICM wind direction and the plane of the disc. 

4.2. 2D velocity fields of ram-pressure stripped galaxies 

With the help of 2D velocity fields the nature of the interaction 
becomes more accessible than with simple rotation curves from 
long-slit spectroscopy. Quantitative analyses on the kinematic 
axis, the kinematic centre and more detailed analyses such as 
harmonic expansion, which is, for example, used by the kineme- 
try package of Krajnovic et al. (2006), offer additional possi- 
bilities to identify distortions in a VF. We have already applied 
the kinemetry method to tidally disturbed simulated systems in 
Kronberger et al. (2007). 

In Fig. |2] we show such a 2D velocity field of a face-on 
ram-pressure affected galaxy after 250 Myr of ram pressure act- 
ing. (The distortions at this time step are qualitatively equal to 
those present after 100 Myr, but more prominent. Therefore we 
chose this particular time step, without loss of generality). The 
distortions are not as severe as those caused by tidal interac- 
tions, which were presented in Kronberger et al. (2007). An al- 
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Fig. 5. Rotation curve of a model galaxy after 100 Myr of ram 
pressure acting face-on (top panel) and after 500 Myr of ram 
pressure acting edge-on (bottom panel). Both sides of the RC 
are plotted on one side for better judgement of the asymmetry. 
The blue-shifted part of the RC is plotted as blue solid line, while 
the redshifted part is plotted as dashed red line. 



most undisturbed inner part of the velocity field is clearly visi- 
ble, which we indicated with an ellipse. The semi-major axis of 
this ellipse corresponds roughly to the stripping radius (12 kpc) 
of this model system as estimated in Kronberger et al. (2008). 
Beyond this radius the VF appears disturbed. This picture of an 
undisturbed inner part and a severely affected outer part beyond 
the stripping radius also agrees with the rotation curve of the sys- 
tem presented in the previous section (c.f. Fig.|2]). As for the ro- 
tation curve, the VF appearance depends on the inclination. The 
characteristic features become less prominent when the galaxy 
is observed close to edge-on. In Fig. [8] we study the appearance 
of this velocity field when put at intermediate redshift (the same 
investigation as performed in Kronberger et al. (2007) for tidally 
disturbed velocity fields). 

The VF appears smooth and widely undisturbed at z = 0.4, 
as the irregularities are smeared out, mainly due to seeing (there- 
fore the use of integral field spectrographs in combination with 
adaptive optics is desirable). The radial profiles of the kinemetric 
properties, which we calculated using the kinemetry programme 
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Fig. 6. Asymmetry measure as a function of the angle between 
the disc and the ICM wind direction (0° is edge-on and 90° is 
face-on ram pressure). The asymmetry measure was averaged 
over several viewing angles. The error bars indicate the standard 
deviation of the mean. 
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Fig. 7. 2D velocity field of a face-on ram-pressure affected 
galaxy after 250 Myr of ram pressure acting. The almost undis- 
turbed inner part of the velocity field is highlighted by an ellipse, 
whose semi-major axis roughly corresponds to the stripping ra- 
dius (~ 12 kpc as estimated in Kronberger et al. 2008). The centre 
of the stellar disc is indicated by the black cross. The black arrow 
in the top right corner indicates the projected ICM wind direc- 
tion. Note that the ICM wind acts perpendicular to the plane of 
the disc. 



of Krajnovic et al. (2006), are presented in Fig.|9]for the model 
galaxy after 250 Myr of face-on ram pressure acting. The posi- 
tion angle F and the flattening q of the best fitting ellipses as well 
as the first and the fifth order Fourier terms ki and k;, are plotted 
as a function of radius for redshift z = 0.1 (top) and z - 0.4 (bot- 
tom). In the radial behaviour of these quantities the undisturbed 
inner part of the velocity field and the distorted outer part are 
visible. At a radius of ~ 5.5 " (at z=0.1, corresponding to ~ 10 
kpc) a clear increase in the ratio of k<, and k[ can be seen. This 
ratio represents complex, kinematically separate components in 
the velocity field. At the same radius the flattening of the el- 
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35 kpc 



Fig. 8. The 2D velocity field of a face-on ram-pressure affected 
galaxy after 250 Myr of ram pressure acting as seen at redhift 
z-0.1 (top panel) and z-0.4 (lower panel), respectively. The cen- 
tre of the stellar disc is indicated by the black cross. The black 
arrow in the top right corner of each panel indicates the projected 
ICM wind direction. Note that the ICM wind acts perpendicular 
to the plane of the disc. 



lipses and the first order Fourier term ki also show significant 
variations, which are, however, lower than in the case of tidally 
induced distortions (c.f. corresponding figures in Kronberger et 
al. 2007). At a redshift of z = 0.4, the variations of the radial 
kinemetric properties get less prominent. A clear identification 
of ram pressure affected galaxies is especially difficult as at such 
redshifts often only the inner part of the velocity field is observ- 
able. In the outer regions of the disc some irregularities in the 
velocity field remain visible also at intermediate redshift. They 
are observable as increase in the ratio of ks and ki . 

In Fig. [To] we present the 2D velocity field of an edge-on 
ram-pressure affected galaxy after 400 Myr of ram pressure act- 
ing as seen at redhift z-0. 1 (top panel) and z=0.4 (lower panel), 
respectively. The mismatch between kinematic centre and the 
centre of the stellar disc is the strongest effect of edge-on ram 
pressure. Already by visual assessment of the VF presented in 
Fig. [TO]several distortions are visible also in the inner part of the 
VF. These iiTegularities are, however, completely smeared out at 
redshift z - 0.4. Also the radial profiles of the kinemetric prop- 
erties, presented in Fig.[TT|reflect this behaviour. At z = 0. 1 vari- 
ations of the flattening parameter q and of the ratio of k^ and ki 
are already visible at smaller radii. The amplitude of these vari- 
ations is however small. To study the spatial distribution of the 
cold gas seems therefore to be still the best way to identify ram 




t i - . * ' 1 

0,5 1 1.5 2.0 2.5 

R [arcsec] 



Fig. 9. Radial profiles of the kinemetric properties, calculated 
using the kinemetry programme for the model galaxy after 250 
Myr of face-on ram pressure acting, for redshift z = 0.1 (top) 
and z - 0.4 (bottom). The position angle F and the flattening q 
of the best fitting ellipses as well as the first and the fifth order 
Fourier terms ki and ^5 are plotted as a function of radius. 



pressure acting in a galaxy. The influence on the velocity field is 
relatively low, as the rotational velocity depends mostly on the 
gradients of the gravitational potential and the gas accounts only 
for a small mass fraction. 

5. Discussion 

Ram-pressure stripping has been modelled numerically using 
different techniques. All these approaches use some simplifi- 
cations due to the complexity of the multi-phase, multi-scale 
gas physics present in the inter-stellar medium of galaxies. One 
method is to treat the ISM hydrodynamically either with an 
Eulerian code (e.g. Roediger & Briiggen 2007) or using SPH 
(e.g. Jachym et al. 2007). SPH has the advantage that it can 
be easily combined with collisionless N-body simulations for 
the dark matter and stellar component. In the setup used for 
this work the collisionless component and the gas are treated 
fully self-consistently. An alternative method is to model the 
ISM with particles which can have inelastic collisions, so-called 
'sticky particles (e.g. VoUmer et al. 2006). The ICM ram pres- 
sure is then modelled as additional acceleration on particles. 
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35 kpc 



Fig. 10. 2D velocity field of an edge-on ram-pressure affected 
galaxy after 400 Myr of ram pressure acting as seen at red- 
hift z=0. 1 (top panel) and z=0.4 (lower panel), respectively. The 
black arrow in the top left of each panel indicates the ICM wind 
direction and the cross and the circle indicate the kinematic and 
stellar disc centre, respectively. 



IF 350 




R [arcsec] 

Fig. 11. Radial profiles of the kinemetric properties, calculated 
using the kinemetry programme for the model galaxy after 400 
Myr of edge-on ram pressure acting, for redshift z = 0.1 (top) 
and z - 0.4 (bottom). The position angle F and the flattening q 
of the best fitting ellipses as well as the first and the fifth order 
Fourier terms ki and are plotted as a function of radius. 



which are not shielded by others. All these methods have their 
advantages and restrictions. If the ISM is considered as contin- 
uous fluid the hydrodynamic prescription yields a better model 
than sticky particles. The latter method is more suitable for non- 
continuous systems (e.g. VoUmer et al. 2006). The SPH method 
has intrinsic shortcomings in the resolution of instabilities and 
turbulence (e.g. Agertz et al. 2007). In systems affected by ram- 
pressure stripping shear flows are present, which can lead to the 
development of Kelvin-Helmholtz instabilities (e.g. Roediger 
et al. 2006). These are not resolved in our setup. The Kelvin- 
Helmholtz instabilities affect mostly the stripped material in the 
wake, where the shear flow is strong and the gravitational force 
is weak (Roediger et al. 2006). The effect for the bulk motion 
of the gas in the plane of the disc is probably small. A similar 
study as the one presented in the present work should be car- 
ried out with a Eulerian code in a future work to investigate the 
importance of hydrodynamic instabilities in this issue. The huge 
advantage of our setup is that we treat gas, dark matter and the 
stellar component fully self-consistently. 

Some of the results presented in this work were also found 
in previous studies. Roediger & Briiggen (2007), for example, 
found that the gas distribution of ram-pressure affected galax- 
ies is more asymmetric for edge-on ram pressure. Consistently 



we find this trend also in the velocity field. The mismatch be- 
tween the kinematic and the stellar disc centre is caused by two 
effects. Ram pressure affects the whole disc. The gas gets dis- 
located where the restoring gravitational force is too low. This 
ram-pressure pushing (Roediger & Briiggen 2007) leads also 
to a central density enhancement and consequently to an in- 
creased central star formation (Kronberger et al. 2008; VoUmer 
et al. 2001). Within our model this enhanced star formation is 
self-consistently treated and induces an increased galactic wind. 
These feedback mechanisms, which do affect the gas kinematics 
due to modified density and pressure gradients, were neglected 
in previous studies and are one of the strengths of the present 
work. With an offset of ^2 kpc, the effect is, however, relatively 
small. 

In a series of papers Vollmer et al. use the velocity field infor- 
mation to fit models of ram-pressure affected galaxies to HI ob- 
servations of nearby cluster galaxies (e.g. Vollmer et al., 2008). 
They identify characteristic features in the VFs of the observed 
systems and study in which configuration the model can repro- 
duce these features. In Vollmer et al. (2008) the velocity field of 
an almost edge-on stripped galaxy is shown. The VF is quite reg- 
ular as the galaxy is in an early phase of stripping. Thus a direct 
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comparison to our models is difficult. Even for edge-on stripping 
the VF distortions are moderate in our models. This is also true 
for the simulations of Vollmer et al. (2008) (compare e.g. their 
Fig. 20). Also in Vollmer et al. (2006) the VF appears rather reg- 
ular. The VF of NGC 4569 presented in Vollmer et al. (2004) 
shows an asymmetric velocity field with an angle between or- 
bital and disc plane of 35°. This is consistent with our finding 
that the velocity field becomes asymmetric for angles <45°. 

More detailed comparisons with observed optical velocity 
fields and rotation curves from the sample of Ziegler et al. (2007) 
will be subject of an upcoming work. 



6. Summary and conclusions 

We have investigated 2D velocity fields and rotation curves 
of the gas in ram-pressure affected spiral galaxies using N- 
body/SPH simulations. 

- We found that the appearance of distortions of the gaseous 
disc due to ram-pressure stripping depends on the direction 
of the acting ram-pressure. In our investigation we studied 
the two extreme cases of face-on and edge-on ram pressure. 
In both cases, the velocity fields of the ram-pressure affected 
galaxies show characteristic features. 

- In general, the distortions of the gas velocity field are sig- 
nificantly lower than in the case of distortions induced by 
tidal interactions. The reason is that the rotational velocity 
depends mostly on the gravitational potential and the gas ac- 
counts only for a small mass fraction. 

- In the case of face-on acting ram pressure distortions appear 
mainly in the outer parts of the galaxy. The separation be- 
tween undisturbed inner and symmetrically disturbed outer 
part is visible in both the rotation curve and the 2D velocity 
field. 

- In the case of ram pressure acting edge-on we find a sig- 
nificant mismatch between kinematic centre and the centre 
of the stellar disc. The rotation curve is asymmetric and the 
outer parts of the RC are declining. The 2D velocity field 
shows signatures of the interaction also in the inner part of 
the disc. 

- At an angle of 45° between the ICM wind direction and 
the disc, the velocity field asymmetry increases significantly 
with respect to smaller angles. The scatter is, however, in- 
creasing towards the edge-on wind too. Therefore it seems 
difficult to distinguish further between different ICM wind 
directions purely with the help of the VF asymmetry. 

- The coUisionless stellar disc is not affected by ram pressure 
and hence shows no disturbed kinematics. The characteris- 
tic features in the gas kinematics presented here, together 
with an absence of distortions in the stellar kinematics are 
therefore a strong indication for the presence of ram-pressure 
stripping acting on a cluster galaxy. 
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